Electron coherent and incoherent pairing instabilities in inhomogeneous 
bipartite and nonbipartite nanoclusters 
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Exact calculations of collective excitations and charge/spin (pseudo)gaps in an ensemble of bipar- 
tite and nonbipartite clusters yield level crossing degeneracies, spin-charge separation, condensation 
and recombination of electron charge and spin, driven by interaction strength, inter-site couplings 
and temperature. Near crossing degeneracies, the electron configurations of the lowest energies 
control the physics of electronic pairing, phase separation and magnetic transitions. Rigorous con- 
ditions are found for the smooth and dramatic phase transitions with competing stable and unstable 
inhomogeneities. Condensation of electron charge and spin degrees at various temperatures offers 
a new mechanism of pairing and a possible route to superconductivity in inhomogeneous systems, 
different from the BCS scenario. Small bipartite and frustrated clusters exhibit charge and spin 
inhomogeneities in many respects typical for nano and heterostructured materials. The calculated 
phase diagrams in various geometries may be linked to atomic scale experiments in high Tc cuprates, 
manganites and other concentrated transition metal oxides. 
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I. INTRODUCTION 

Strongly correlated electrons in cuprates, manganites 
and other transition metal oxides exhibit high Tc su- 
perconductivity, magnetism and ferroelectricity accom- 
panied by spatial inhomogeneities at the nanoscale level 
Ul) S B 0) S B 0) Hi ■ Over the past few years there is an 
increase in interest to electron instabilities in nanoclus- 
ters, assembled clusters of correlated materials in vari- 
ous topologies for synthesizing new nanomaterials with 
unique electronic and magnetic properties 
viously, there is a clear need for an accurate analysis 
of electron correlations, fluctuations and instabilities in 
nanoclusters and large complex systems with competing 
phases. The closed form solution, existing in the Bethe 
ansatz ground state [ll|, is difficult to analyze at finite 
temperatures T > without having to resort to various 
approximations. Perturbation theory is usually inade- 
quate while numerical methods have serious limitations, 
such as in the Quantum Monte Carlo method with its no- 
torious sign problem where the resulting approximations 
often lead to some controversy. On the contrary, exact 
calculations in small clusters [l^, [H, [fl E^j [3 S^^^ 
appealing alternative for the detection of possible phase 
separations and spatial inhomogeneities especially at fi- 
nite temperatures. As far as the authors are aware, an 
exact analysis of level crossing instabilities (degeneracies) 
in canonical ground state eigenvalues and correspond- 
ing competing average energies at finite temperature for 
a general on-site interaction U and electron concentra- 



tions have not been attempted in small or moderate size 
clusters |l7[. Exact computations of electron instabili- 



ties in various cluster geometries at the nanoscale level 
can be vital to the understanding of the role of thermal 
and quantum fluctuations for large pairing gaps and a 
transition temperature Tc in the correlated nanoclusters, 
nanomaterials and corresponding "large" inhomogeneous 
systems [11 El HO, HH, M, MMM- 

Although our approach for "large" systems is only ap- 
proximate, this class of clusters in an ensemble displays 
a common behavior which we believe is generic for large 
thermodynamic systems. Our results for typical bipar- 
tite and frustrated (nonbipartite) cluster geometries have 
successfully mapped out scenarios where many body lo- 
cal effects are sufficient to describe spin-charge separa- 
tion and pairing pseudogaps at the nanoscale level. Spa- 
tial microscopic inhomogeneities have been observed in a 
number of scanning tunelling microscopy (STM) probes 
in doped high-Tc superconductors (HTSCs). There is 
growing evidence suggesting that inhomogeneities at the 
nanoscale level, in the so-called stripes surrounded by 
essentially neutral correlated MH-like antiferromagnetic 
insulators [1^, [l^], play a defining role for the electron 
pairing and the orig in of superconductivity at the atomic 
scale in HTSCs [28l . [29| . Besides the existence of charge 
pairing, the inhomogeneities of possible electronic nature 
can exist in a form of spatially separated magnetic phases 
in cuprates and manganites under doping [30| • The mag- 
netic inhomogeneities seen in other transition metal ox- 
ides at the nanoscale level, widely discussed in the lit- 
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erature [3l|, |32, l33|, |3J, |35|, can be crucial for the spin 
pairing instabilities, origin of ferromagnetism and ferro- 
electricity in the spin and charge subsystems [s^, [s^l ■ A 
phase separation of the ferromagnetic clusters embedded 
in an insulating matrix is believed to be essential to the 
colossal magnetoresistance (CMR) in manganese oxides. 
At sufRciently low temperatures, the spin redistribution 
in an ensemble of clusters can produce inhomogeneities 
in the ground state and at finite temperatures [2J| . The 
non-monotonous behavior of the chemical potential ver- 
sus electron concentration found in generalized self con- 
sistent approximation [sll also suggests possible electron 
instabilities and inhomogeneities near half filling. From 
this perspective, exact studies at T > of electron charge 
and spin instabilities at various U > 0, inter-site cou- 
plings and various cluster topologies can give important 
clues for understanding of charge/spin inhomogeneities 
and local deformations for the mechanism of pairings 
and magnetism in "large" concentrated systems when- 
ever correlations are local. 

It is a generally believed that a strong on-site Coulomb 
interaction supports ferromagnetism and is detrimental 
for the electron pairing and superconductivity in clusters 
and "large" concentrated systems [s^. Our exact stud- 
ies of gaps and pseudogaps in finite-size systems have 
uncovered some important answers related to spin-charge 
separation, pairing and thermal condensation of the elec- 
tron charge and spin. Despite this, there is still a vast 
amount of uncertainties that need to be unravelled: (i) 
What are the conditions for the electron phase separa- 
tion instabilities and spin/charge inhomogeneities? (ii) 
What is the role of inhomogeneities and are these spa- 
tial spin/charge inhomogeneities crucial for the pairing 
mechanisms in these compounds? (iii) When treated ex- 
actly, what essential features can the Hubbard clusters 
capture that share similar properties with the "large" 
concentrated transition metal oxides? 

A redistribution of excess electron/hole inhomo- 
geneities or spin up/spin down domains in an ensemble 
of tetrahedrons for all ?7 > depends on the sign of the 
hopping term [23|. Here we show that in the distorted 
square pyramids in the perovskite structures, the inter- 
site coupling c between the apex site with the base can be 
beneficial or detrimental for the electron pairing or fer- 
romagnetism. An unstable "saturated ferromagnetism" , 
existing in frustrated lattices at low temperatures and 
large U for a particular sign of hopping {t > 0) [11] , im- 
plies either antiferromagnetism, unsaturated ferromag- 
netism, or electron coherent pairing for charge and spin 
pairing (pseudo)gaps. Here it is argued that for one hole 
off half filling electrons undergo separate thermal con- 
densation of the charge and spin degrees (independent of 
cluster topology); the system may be divided into two 
coexisting and dynamically bound bosonic subsystems, 
where two types of individual bosonic pairs, made up of 
double electron charges and oppositely oriented (antipar- 
allel) spins, can fluctuate. We shall see that the phase 
diagram, under some circumstances, is mostly controlled 



by the changes in the cluster geometry (topology). 



II. MODEL AND FORMALISM 



It is possible to assume that the electron pairing and 
magnetic instabilities of the purely electronic nature is 
described by a local Coulomb interaction J7 in a single 
band Hubbard model 



H 



E 



(1) 



The sign of the hopping amplitude t between the nearest 
neighbor sites in ([1]) leads to essential changes of elec- 
tronic structure. In nonbipartite clusters, such as tetra- 
hedron, we considert = ±1. In addition, for the distorted 
pyramid we take the coupling parameter between the api- 
cal site and the atoms in the base equal to ci, with c < 1. 
Our studies of the quantum and thermal fluctuations of 
electrons in finite clusters are based on exact diagonal- 
ization, analytical and numerical calculations of energy 
levels and expressions for the canonical and grand canoni- 
cal partition functions in various cluster geometries. The 
exact grand canonical potential for the interacting 
electrons {U) in an external magnetic field {h) is 



fir 



-Tln^:, 
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where N and are the number of particles and the pro- 
jection of the spin in the n-th quantum state. The first 
and second order responses of the charge and spin degrees 
due to the changes in the chemical potential /i (doping) or 
an applied magnetic field are calculated without taking 
the thermodynamic limit. The competing energy states, 
in conjunction with the canonical and the grand canoni- 
cal ensemble, yield valuable insight into electron instabil- 
ities in the real nanoclusters and nanomaterials with the 
correlated electrons. The introduced formalism allows 
us to describe the smooth and sharp phase transitions 
with competing stable and unstable inhomogeneities in 
the canonical and grand canonical ensembles. 

Below in Sec. IIIII we provide a detailed description of 
the general methodology: we define the criteria for the 
charge and the spin pairing instabilities in the canonical 
and grand canonical ensembles; formulate the conditions 
for existence of quantum critical points, coherent pairings 
and spontaneous transitions in the ground state and cor- 
responding critical temperatures of crossovers for various 
phases and boundaries in the phase diagrams discussed 
in Sees. IIV] and IVl 



III. GENERAL METHODOLOGY 
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A. Canonical charge and spin gaps 

To facilitate the comparison with the frustrated clus- 
ters, we summarize here the main results in the ground 
state and at the finite temperatures for bipartite and non- 
bipartite clusters obtained earher in Refs. p^. [20l. I21I |2^. 
[23, m, . The degrees of freedom for charge and spin, 
electron and spin pairings, temperature crossovers, quan- 
tum critical points, etc. were extracted directly from the 
thermodynamics of these clusters. One can classify the 
charge and spin order parameters as an energy difference 
between the various competing phases by analogy with 
phase transitions in the thermodynamic limit. In the 
ground state, the calculated differences in the canonical 
energy levels between configurations with various num- 
bers of electron charge and spin determine the energy 
gaps for electron charge and spin excitations. Using the 
exact partition function in the canonical ensemble, we 
also analyzed analytical expressions for the average ener- 
gies for various number of electrons N . For given temper- 
ature T and C/, we calculated the energy differences yLt+ = 
E{N+1)-E{N) and = E{N)-E{N-1) for the aver- 
age canonical energies E{N) by adding or subtracting one 
electron (charge) in the cluster for a given spin S. The 
energy difference between the two consecutive excitation 
energies by adding or subtracting electron can serve as a 
natural order parameter in a canonical approach. Then 
the charge gap at finite temperature can be written as 
A%T) = fj,_ = E{N + + E{N - 1) - 2E{N). 

The opening of the gap is a local correlation effect, and 
clearly does not follow from long range order, as exempli- 
fied here. The difference /i+ — /x_ is somewhat similar to 
the difference I~A for a cluster, where / is the ionization 
potential and A the electron affinity. For a single "impu- 
rity" at half filling and T = 0, I — A is equal to U, which 
represents a screened local parameter U in the Hubbard 
model [40*1 (HI). Thus the gap picture is analogous to an 
inter-configuration energy gap for the crossover between 
different many body ground state ionic configurations in 
solids. For example, the charge gap is simply equivalent 
to the energy of the "reaction" between different cluster 
configurations {d) at fixed N 

d^ + d^ ^d^+^ +d^-\ (3) 

i.e., the difference in the canonical energies of ioniza- 
tion and affinity for many body cluster configurations in 
ensemble. However, the configurational change in the 
ensemble of isolated clusters is supposedly due to the 
possible spontaneous fluctuations in the electron num- 
bers and electron redistribution via a charge reservoir. 
The negative spin gap in the canonical ensemble can 
be treated correspondingly. We calculate a spin gap 
as the difference in the average energies between the 
two cluster configurations with various spin S states, 
A'{T) = E{S E{S), for E{S) being respectively 

the average canonical energy in the spin sector at fixed 
N p. 



B. Charge and spin instabilities 

Many phenomena and phase transitions invoked in the 
approximate treatments of "large" concentrated systems 
are seen also in the exact analysis of pairing instabilities 
in the canonical ensemble of the small clusters in thermo- 
dynamic equilibrium [H, ^ HH HI, [13, [13, (H] . As we 
shall see, in some circumstances, small changes of the ex- 
ternal parameters can lead to level crossing instabilities 
in various electron configurations with the formation of 
negative charge and spin gaps. Physically, a positive gap 
manifests the phase stability and smooth crossover, while 
a negative gap describes spontaneous transitions from 
one stationary state to another. Instead of a full phase 
separation at A'^'^ < 0, the local inhomogeneities in the 
clusters can provoke electron redistribution and quantum 
mixing of the various charge and spin configurations. In 
the presence of a negative gap, the many-body ground 
state has an appreciable probability of being found in 
either of these competing configurations. The collective 
particle excitations are also reflected in the fluctuations of 
the pair density in Eq. It is intriguing that these fluc- 
tuations make the pair redistribution across the clusters 
possible even without direct contact between the clus- 
ters. These fluctuations play a crucial role of the pair 
transitions in the absence of electron hopping between 
clusters in Eq. ((T|). Near ground state degeneracies, the 
lowest energy states control the low energy dynamics of 
the electronic and magnetic transitions over a significant 
portion of the phase diagram. 

The possible quantum critical points, phase transitions 
and nonzero temperature crossovers are described using a 
simple cluster approach: we define critical parameters for 
the level crossing degeneracies or quantum critical points 
from the vanishing conditions for the canonical charge 
and spin gaps, i.e., A'^'^{U, c) = 0. The sign of the gap is 
also important in identifying the regions for the electron 
charge and spin instabilities, such as the electron-electron 
A'^ < 0, electron-hole A° > pairings in the charge sec- 
tor or the parallel A" < and opposite A^ > spin 
pairings in the spin sector. The key question here is the 
exact relationship between the canonical charge A'^ gap 
and its corresponding grand canonical spin A* counter- 
part calculated for various bipartite and frustrated clus- 
ter topologies. For charge degrees the negative sign of gap 
implies phase (charge) separation [i.e., segregation) of 
the clusters into hole-rich (charge neutral) and hole-poor 
regions. The quantum mixing of the closely degenerate, 
hole-poor d^~^ and hole-rich d^^^ clusters for one hole 
off half filling, instead of causing global phase separation, 
provides a stable spatial inhomogeneous medium that al- 
lows the pair charge to fluctuate. The inhomogeneities 
favored by the negative gaps are essential for providing 
the spontaneous redistribution of the electron charge or 
spin. The inhomogeneities in the charge redistribution 
for A'^ < and A* — imply static heterostructure 
for different electron conflgurations, close in energy, in 
an unstable ensemble of clusters. These inhomogeneities 
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are consistent with nucleation of the "negative" charge 
gap in cuprates above Tc 0. At low temperatures, the 
dynamic picture for pair fluctuations between different 
electron configurations A'^ < and A" > is possible at 
relatively low temperatures in spatially inhomogeneous 
coherent state (A" = -A^, see also Sec. HV^ . This 
result is consistent with the observation of nonlocal su- 
perconductivity at low excitation energies and at higher 
energies, holes localized in an inhomogeneous "stripe" 
pattern The negative spin gap describes the possible 
parallel spin pair binding instability. This picture im- 
plies spontaneous ferromagnetism and phase (spin) sepa- 
ration into domains in accordance with the Nagaoka theo- 
rem. For the negative gaps, one can introduce the critical 
temperatures (fi) and Tf (/i) versus chemical poten- 
tial for boundaries between various phases derived from 
the condition that the corresponding gaps disappear, i.e., 
A'='"(T,/x) = 0. 



C. Charge and spin susceptibility peaks 

Conventional phase transitions at finite tempera- 
ture are driven by thermal fluctuations. In the grand 
canonical approach using exact analytical expressions 
for the grand canonical potential and partition func- 
tions as expressed in Eg. pi), we have analyzed (in 
Refs. [H ia mi 111 m iMlil) the variation of the 
charge, and spin, density of states or correspond- 
ing charge Xc(m) and spin Xs(/^) susceptibilities. 



Xc 



d(N) 



Xs 



dh 



(4) 



as a function of the chemical potential /i and /i in a wide 
range of temperatures. In a grand canonical approach 
the energy difference between the two consecutive sus- 
ceptibility peaks in terms of /i and h at finite tempera- 
tures can serve as a natural order parameters for charge 
and spin degrees respectively. This energy difference for 
density of states in /i space determines the charge gap 
in canonical approach. We find (opposite) spin pairing 
gap by calculating the minimal magnetic field necessary 
to overturn the spin. In the grand canonical method 
we define the gap as a magnetic field at which the dis- 
tance between the subsequent spin susceptibility peaks 
in fi space vanishes. Using the maxima of zero mag- 
netic field susceptibility, ^\h->o, we also calculated the 
boundary curve for the onset of the spin gap for vari- 
ous in infinitesimal ft, — > above . To distinguish 
this from the canonical and grand canonical gaps at finite 
temperatures we call it pseudogap. The opening of such 
distinct and separated (pseudo)gap regions for the spin 
and charge degrees at various fillings in /i space is indica- 
tive of the corresponding spin-charge separation. The 
crossover temperatures and phase boundaries for vari- 
ous transitions can be found by monitoring maxima and 
minima in charge and spin susceptibilities. We define 
the critical temperatures Tc and T* in equilibrium as the 
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FIG. 1: Thermal occupation probabilities versus fj. in the 
grand canonical ensemble of 4-site clusters in the vicinity of 
quantum critical point fip for (A'^) « 3 at [/ = 4.0. The phase 
separation at relatively low temperatures below T < 0.0075 
manifests a significant suppression of (A*') — 3 clusters close to 
optimal doping fip = 6.557. In this area (A'^) — 2 and (A^) — 4 
clusters share equal weight probabilities, while at higher tem- 
peratures (A'^) « 3 also becomes thermodynamically stable. 
In equilibrium, the grand canonical value fi p a,t optimal dop- 
ing at r = reproduces the result i-ip = + At-)/2 for the 
canonical approach. 



temperature at which the distances between the charge or 
spin susceptibility peaks vanish and corresponding pseu- 
dogaps disappear (see Sec. IV Bp . Notice that according 
to the given definition, the energy pseudogaps obtained 
in the grand canonical method are positive which is a key 
difference from the canonical gaps. 



D. Charge and spin inhomogeneities 

The developed grand canonical approach can be ap- 
plied to understand of the electron fluctuations and the 
spatial inhomogeneities to model the behavior of the con- 
centrated systems in bipartite and frustrated structures. 
An ensemble of bipartite clusters at small and moderate 
U exhibits typical inhomogeneous behavior in its charge 
distribution. A normalized probability wjv for the elec- 
tron distribution in grand canonical ensemble as a func- 
tion of temperature T for various electron numbers N is 
the following 



71, N 



(5) 



The calculated probabilities of electrons in competing 
conflgurations are shown in Fig. [T] for the 4-site cluster 
at t/ = 4. At low temperatures and electron concentra- 
tion close to /Ltp, the clusters with (TV) = 2 and {N) = 4 
have equal probabilities, W2 = W4 « 0.5. In some circum- 
stances electron conflgurations in equilibrium can have 
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close energies for the clusters in contact with a particle 
reservoir. This picture shows a mixture of ungapped and 
partially gapped states. As temperature increases, the 
probability UJ3 of (N) = 3 clusters with unpaired spin 
gradually increases, while the probability of finding spin 
paired, hole-rich and hole-poor clusters decreases. 

Qualitatively, the formation of inhomogeneous electron 
distribution or "stripe" picture can be understood from 
simple energy considerations (see Sec|n|. For a fixed av- 
erage number of electrons, the charge and spin on each 
separate cluster in the ensemble can fluctuate. The two 
configurations close in energy are nearly degenerate, and, 
as temperature increases, it is energetically favorable to 
have some clusters with d^~^ and another with 
instead of having clusters with electrons. These re- 
sults, that depend on the cluster geometries, parameter 
U as well as on the sign of t, can be directly applied 
to nano and heterostructured materials, which usually 
contain many independent clusters, weakly interacting 
with one another with the possibility of having inhomo- 
geneities for a different number of electrons per cluster. 
At half filling, the antiferromagnetic state has the lowest 
energy per electron. Therefore, the energy can be min- 
imized upon small doping by segregration of holes into 
charged clusters with different number of electrons. The 
embedded antiferromagnetic background with opposite 
spin pairing provides a spin rigidity (unperturbed) me- 
dia that allows inhomogeneities to optimize the coherent 
pair fluctuations across the clusters [l^l- The mixture 
of the closely degenerate ferromagnetic domains can also 
lead to the stable spatial magnetic inhomogeneities for 
spin fluctuations. Interestingly, the quantum and ther- 
mal fluctuations in the canonical and grand canonical en- 
sembles display "checkerboard" patterns nanophase 
inhomogeneities [26| and nucleation of pseudogaps driven 
by temperature seen recently in nanometer and atomic 
scale measurements ^ in HTSCs above Tf [1,0]. Mi- 
croscopic spatial inhomogeneities and incoherent pairing 
pseudogaps in nanophases measured by scanning tun- 
nelling microscopy (STM) correlate remarkably with our 
predictions using small 4-site and 2x4 nanoclusters [13] ■ 



E. Coherent charge and spin pairings 

The behavior of such clusters near crossing degenera- 
cies in a quantum coherent phase with minimal spin at 
low temperatures is somewhat similar to the conventional 
BCS superconductivity (see Sec. IIV A]) . We found that 
at rather low temperatures the calculated positive pseu- 
dospin gap A* in the grand canonical method can have 
equal amplitude with a negative charge gap A'^ derived 
in canonical method. A* = ]A'^]. Such behavior is similar 
to the existence of a single gap in the conventional BCS 
state. We call such an opposite spin (singlet) couphng 
and electron c harg e pairing as a spin coherent electron 
pairing in Ref. [2J]. However, unlike to the BCS theory, 
the charge gap differs significantly from the spin pseu- 




FIG. 2: Charge A'^ and spin A** gaps versus U in an ensemble 
of tetrahedrons at t = 1, (N) « 3 and T = 0.001. Negative 
charge gap A'^ < implies charge phase separation, while 
the positive, opposite spin pairing gap of equal amplitude 
A" = —A'^ describes Bose condensation of electrons similar 
to BCS-Iike coherent pairing with a single (unique) energy 
gap. This coherent state is analogous to Phase A in 4-site 
clusters [2^ . The spin gap has been calculated using the 
grand canonical approach. 




1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 
U 



FIG. 3: The charge A'' > and parallel (triplet) spin A^ < 
gaps versus U in an ensemble of tetrahedrons at t = — 1, 
(A^) — 3 and T = 0.001. The positive charge gap for all 
U describes MH-Iike insulating behavior analogous to Phase 
C in 4-site clusters [2^. The negative spin gap (A^ < 0), 
coexisting with the charge pairing gap displays S = ^ Na- 



gaoka saturated ferromagnetism at all f/ > (see Sec. IIV B|l . 
The charge gap has been calculated using grand canonical 
approach. 
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dogap as temperature increases above . For example, 
the vanishing of double peak structure in zero spin sus- 
ceptibility gives a critical temperature , at which the 
spin pseudogap disappear. The canonical charge gap dis- 
appears at higher temperatures, i.e., A'^{T^) = 0. The 
BCS-like coherent behavior and possible superconduc- 
tivity with condensation of opposite spin pairs occur at 
rather low temperatures (see Sec. IIIICp . while electron 
charge pairing can be established at relatively high tem- 
peratures, < . The positive spin gap calculated 
in grand canonical approach implies homogeneous elec- 
tron (opposite) spin spatial distribution below Tf. This 
picture is consistent with the spatially homogeneous spin 
pseudogap that opens below Tc = Tf in doping depen- 
dent STM measurements of Bi2Sr2Cu06+a; Q- We also 
find a close analogy for the coherent electron pairing in 
clusters with real space singlet pairs in resonance valence 
bond states or local inter confi gur ation fluctuations in 
mixed valence states [HI, [H, \MMi, El Ell- 

IV. GROUND STATE PROPERTIES 



A. Bipartite clusters 

Exact calculations for charge and spin gaps in small 
clusters in various geometries are important for under- 
standing the electron ground state behavior in bipartite 
and nonbipartite (frustrated) systems. Below we sum- 
marize the results for electron instabilities and phases 
obtained earlier (see Fig. 1 in Ref. [23|) for square and 
other bipartite clusters with one hole off half filling at 
infinitesimal T — > 0. The vanishing of gaps at quan- 
tum critical points, Uc — 4.584 and Up — 18.583, in- 
dicates energy level crossings and electron instabilities 
in 4-site clusters for charge and spin, respectively. The 
charge A'^ and spin A* gaps versus U in an ensemble 
of square clusters at (N) « 3 exhibit at infinitesimal 
T — > the following phases; Phase A: Charge and spin 
pairing gaps of equal amplitude A^ = A^ = — A^ at 
U < Uc describe Bose condensation of electrons similar 
to BCS-like coherent pairing with a single energy gap; 
Phase B: Mott-Hubbard like insulator with A"^ > and 
gapless S — ^ excitations at Uc < U < Up describes 
a spin liquid behavior; Phase C: Parallel (triplet) spin 
pairing (A* < 0) displays 5" = § the saturated ferro- 
magnetism at [/ > C/p in Mott-Hubbard insulator for a 
positive charge gap, A"^ > 0. Notice, that incoherent op- 
posite spin pairing jA''! 7^ A"^, different from the charge 
pairing at U < Uc, suggests spin-charge separation for 
spin and charge degrees at U > Up- 

Square clusters at weak and strong couplings share 
common important features with 2x4 ladders and other 
bipartite clusters fs^]. Negative gaps describe possible 
hole binding or parallel spin pairing instabilities. For 



charge gap vs coupling to vertex 




FIG. 4: Charge gap A'^ versus coupling c between the apex 
atom and the four base atoms in deformed square pyramid 
{t = 1) for one hole of half filling, (iV) = 4, [/ = 3 and 
T — 0.01. Charge and spin pairing gaps of equal amplitude 
= A^ = — A'^ at c < 0.35 imply coherent pairing, while 
A'^ > and A" < at c > 0.35 correspond to a ferromagnetic 
insulator for S = ^. 

charge degrees at weak coupling, this gives an indica- 
tion of phase separation (i.e., segregation) into hole-rich 
(charge neutral) and hole-poor clusters. In contrast, at 
strong coupling the negative spin pairing gap for par- 
allel spins and positive charge gap reveal ferromagnetic 
instability in accordance with the Nagaoka theorem. In 
large bipartite clusters at intermediate U, electrons be- 
have differently from square clusters. For example, in 
2x4 ladders we found an oscillatory behavior of charge 
gap as a function of U [l^ . The vanishing of the charge 
gaps, manifesting the multiple level crossing degenera- 
cies and electronic instabilities in charge and spin sectors, 
is indicative of possible electron instabilities in bipartite 
clusters at moderate U. 



B. Tetrahedrons 

For comparison with small bipartite clusters in 
Sec. lIV Al we consider here a minimal four site nonbipar- 
tite structure. A tetrahedron has a topology equivalent 
to that of a square with the next nearest neighbor cou- 
pling (t' = t) and may be regarded as a primitive unit 
of typical frustrated system. Nonbipartite systems, with- 
out electron-hole symmetry, exhibit a pairing instability 
that depends on the sign of t. Notice that sign of t also 
leads to essential changes in the electronic structure. The 
tetrahedral clusters show pairing instabilities for charge 
degrees at t = 1 and spin degrees at i = — 1 that maxi- 
mizes the amplitudes of negative charge A'^ < and spin 
A*' < gaps and corresponding condensation tempera- 
tures, T^{^) and Tf{^) The negative gap in the 
canonical approach displays electron pairing A^ = |A'^| 
instability for all U. Fig.[3]illustrates the charge and spin 
gaps at small and moderate U. The negative charge gap 
in Fig. [2] is indicative of the inhomogeneous charge re- 
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distribution and phase separation of electron charge into 
hole-rich (charged) and hole-poor (neutral) cluster con- 
figurations [2^. The phase diagram for t = 1 is similar 
to the Phase A in Sec. IIV Al but applied for all U values. 
In contrast, the positive spin gap in the grand canonical 
approach A** > corresponds to uniform opposite spin 
distribution in Fig. O This BCS-like picture for charge 
and spin gaps of equal amplitude A"* = A^ = — A'^ at 
(iV) « 3, in analogy with the square clusters, will be 
called coherent pairing (CP) (2^ . In equilibrium, the spin 
singlet background {xs > 0) stabilizes phase separation 
of paired electron charge in a quantum CP phase. Fig. [2] 
illustrates the charge A'^ and spin A* gaps in tetrahedral 
clusters at {N) « 3, T ^ 0. The unique gap, A* = A^ 
at T = 0, in Fig. [5] is consistent with the existence of 
a single quasiparticle energy gap in the BCS theory for 
U < Q Positive spin gap for all U provides pair 

rigidity in response to a magnetic field and temperature 
(see Sec. IV Bl) . Notice that the coherent pairing exists 
also at large U where Nagaoka theorem for nonbipartite 
clusters with specific sign of t can be applied. The stabil- 
ity of minimal spin 5 = (singlet) state in tetrahedron 
at t = 1 is consistent with that of non maximum (unsat- 
urated) spin in Nagaoka problem. Thus our result shows 
that Nagaoka instability toward spin flip at large U in 
frustrated lattices with t = 1 can be associated with the 
BCS-like coherent pairing applied for general U . 

The negative spin gap A^ in Fig.[3]is shown for canoni- 
cal energy differences between S — ^ and S = \ configu- 
rations. Correspondingly, the positive charge gap A^ > 
in a stable MH-like state is derived using grand canonical 
energies 20j. As in bipartite square clusters, the grand 
canonical positive charge gap A'^ > is different (inco- 
herent) from the parallel spin pairing gap. A* < 0. Thus 
the phase diagram for t = —1 with A"^ ^ —A* is similar 
to the Phase C in Sec. IIV Al but applied for all U val- 
ues in the phase diagram. The negative spin gap for all 
couplings implies parallel spin pairing and Nagaoka-like 
saturated ferromagnetism with maximum spin in the en- 
tire range of U . Spin-charge separation is considered to 
be one of the key properties of the correlated electrons 
that distinguishes t — —1 from t — 1. Such behavior 
at t = — 1 is accompanied by spin-charge separation and 
formation of the rnagnetic (spatial) inhomogeneities or 
domain structures |20| in a wide range of parameters. 



C. Square pyramids 

From the early days of high-temperature supercon- 
ductivity, the idea of a possible role of apical sites in 
p-type superconductors has been controversial. The oxy- 
gen atom position at the apex of pyramidal crystalline 
structure can be altered through the addition of impu- 
rities and can be relocated to a lower or sideways posi- 
tion, thus changing the electron interactions or coupling 
strength c between apex and the planar atoms. There 
is no significant influence of localized electron charge of 



apical site on electron pairing and possible superconduc- 
tivity in CUO2 planes in Bi2Sr2CaCu208+5. When excess 
apex does not exist, i.e., 5 = 0, this system is an insu- 
lator. However, when excess apex oxygen is introduced, 
hole carriers are supplied into Cu02 planes and the ma- 
terial shows superconductivity [48j . Here we try to draw 
a closer connection to HTSCs perovskites and consider 
an ensemble of square pyramids of octahedral structure. 

Fig.|4]shows the charge gap at fixed J7 = 3 and {N) 4 
under the variation of the coupling term c between the 
plane and the apex atoms. This picture gives surprisingly 
plausible evidence for understanding the detrimental role 
of excess electron on charge pairing for possible distor- 
tions of pyramidal crystalline structure in perovskites. In 
Fig. m the strong distortion of the pyramid structure for 
c = (with reduced coordination number) reproduces 
a charge pairing gap in planar square geometries. At 
{N) « 4, the electron is localized and there is no charge 
transfer from apex atom in an ensemble of pyramid clus- 
ters at c = 0. The negative charge gap, identical to 
the spin gap, exists only for c < cq, where cq = 0.35 
is a quantum critical point for level crossing degener- 
acy. Calculated electron distribution, as a function of 
c, shows that electron charge residing on the apical site 
does not contribute to the pairing whenever c is less than 
Cq. The coupling in the pyramid structure at c < cq for 
{N) « 4 leads to charge pairing instability with negative 
charge and positive spin gaps of equal amplitude as seen 
in square clusters at (iV) « 3 in Sec. IIV Al In contrast, at 
c > Co, the induced charge gap driven by c change leads 
to electron hole pairing and a transition into insulating 
Mott-Hubbard (MH) behavior with A'^ > 0. The apex 
atom, coupled to square-planar geometry, have shown to 
have a detrimental affect on the negative charge and posi- 
tive spin gaps, which are favorable to forming a Bose con- 
densate in the region of instability. We found a coherent 
pairing in the phase diagram with one hole off half filling 
also in the ensemble of octahedron clusters (perovskite 
systems) in Ref. There is also shown that octahe- 

dron threaded by magnetic flux in hole-rich regions can 
get trapped in stable minima at half integral units of the 
magnetic quantum flux. Such approach can be applied to 
understand the detrimental effect of the transverse mag- 
netic field on electron charge and opposite spin pairings 
for possible superconductivity in HTSCs in planar face 
centered square (fcs) geometry [23 |. 

V. PHASE DIAGRAM 



A. Tetrahedrons at large U and t = 1 

The charge and spin susceptibility peaks in clusters, 
reminiscent of the singularities in infinite systems, dis- 
play an extremely rich phase diagram at finite tempera- 
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FIG. 5: The T-/i phase diagram of tetrahedrons without 
electron-hole symmetry at optimally doped (A^) « 3 regime 
near fip = 1.998 at U — AQ and t — 1 illustrates the con- 
densation of electron charge and onset of phase separation for 
charge degrees below . The incoherent phase of preformed 
pairs with unpaired opposite spins exists above Ts^ ■ Below 
Tf, the paired spin and charge coexist in a coherent pairing 
phase. The charge and spin susceptibility peaks, denoted by 
T* and Tc, define pseudogap regions calculated in the grand 
canonical ensemble, while phase boundaries fJ,+ {T) and i-i-{T) 
are evaluated in the canonical ensemble. The spin pseudogap 
region exists for Tf < T < T' . Charge and spin peaks rec- 
oncile at T ~ T', while x'^ peak below Tf signifies metallic 
(charge) liquid (see inset for square cluster in Ref. [2^). 

tures. The realization of a high transition temperature, 
Tc, in clusters and bulk systems depends on the interac- 
tion strength U, doping, and the detailed nature of the 
crystal structure (sign and amphtude of t). As exem- 
pUfied here, the critical temperatures for various pairing 
instabilities in frustrated clusters also strongly depend on 
the sign of the hopping {t) term. Fig. [5] for t = 1 illus- 
trates a number of nanophases, defined in Refs. [20l. [2^. 
for the tetrahedron at large U — 40, found earlier in 
tetrahedron and bipartite 2x2 and 2x4 clusters at mod- 
erate /7 = 4 values [l^l- This diagram captures the es- 
sential electron charge and spin pairing instabilities at 
finite temperatures. The curve H+{T) below sig- 
nifies the onset of charge pair condensation. The cal- 
culated susceptibility peaks in Fig. [5] correspond to the 
pseudogap crossover temperature T*. As temperature is 
lowered below T*, a spin pseudogap is opened up first, 
as seen in NMR experiments [22], followed by the grad- 
ual disappearance of the spin excitations, consistent with 
the suppression of low-energy excitations in the HTSCs 
probed by STM and ARPES [l,i,i,[i,l3- In contrast, 
the local charge gap, A'^, evolves smoothly as tempera- 
ture decreases below . The opposite spin CP phase, 
with fully gapped collective excitations, begins to form at 
T <Ts and spin pairing rigidity gradually grows upon 



lowering of the temperature. As temperature decreases 
both charge and spin pseudogaps emerge into one gap 
at zero temperature. Therefore, at sufficiently low tem- 
peratures, this leads to the BCS-like coherent coupling 
of electron charge to bosonic excitations (see Sec. I IV Bp . 
However, the spin gap is more fragile and as temperature 
increases it vanishes at Ts^, while charge pseudogap sur- 
vives until . 

The charge inhomogeneities [l|, in hole-rich and 
charge neutral spinodal regions between fj,^ and fi^ are 
similar to those found in the ensemble of squares and re- 
semble important features seen in the HTSCs. Pairing 
and transfer of holes is a consequence of the existence 
of an inhomogeneous background. In the absence of di- 
rect contact between clusters, the inhomogeneities in the 
grand canonical approach are establishing a transfer of 
paired electrons via this (thermal) bath media. Fig. [S] 
shows the presence of bosonic modes below /Lt+(T) and 
r/' for paired electron charge and opposite spin respec- 
tively. This picture suggests condensation of electron 
charge and spin at various crossover temperatures while 
condensation in the BCS theory occurs at a unique Tc 
value. This result suggests that thermal excitations in 
the exact solution are not quasiparticle-like renormal- 
ized electrons, as in the BCS theory, but collective paired 
charge and coupled opposite spins [2J| • 

The coherent pairing of holes here is a consequence 
of the existence of homogeneous opposite spin pairing 
background, consistent with the STM measurements [29| . 
This led us to conclude that Tf can be relevant to 
the superconducting condensation temperature Tc in the 
HTSCs. In the absence of spin pairing above , the 
pair fluctuations between the two lowest energy states 
becomes incoherent. The temperature driven spin-charge 
separation above Tf resembles an incoherent pairing (IP) 
phase seen in the HTSCs 0, i, S, i, i] . The charged pairs 
without spin rigidity above Tf, instead of becoming su- 
perconducting, coexist in a nonuniform, charge degen- 
erate IP state similar to a ferroelectric phase [2^- The 
unpaired weak moment, induced by a field above Tf, 
agrees with the observation of competing dormant mag- 
netic states in the HTSCs The coinciding x'^ and 
X'^ peaks in the vicinity of critical temperature T' show 
full reconciliation of charge and spin degrees seen in the 
HTSCs above Tc. However, in both channels the charge 
and spin pseudogaps behave differently or independently. 
Indeed, we find that the variation of the spin pairing gap 
with temperature does not cause a change in the charge 
pairing gap. In the absence of electron-hole symmetry in 
the tetrahedrons, the reentrant phenomenon can be ob- 
served at low temperatures ^23] . In Fig. [51 as temperature 
increases near optimal doping ^ < ^p, clusters undergo a 
transition from a CP phase to a MH-like behavior. Notice 
that the charge and spin pairing do not disappear in the 
underdoped regime for fi > iip but are governed predom- 
inantly by the physics of antiferromagnets at half filling. 
In contrast, in the overdoped regime at low temperatures, 
the charge pairing pseudogap gradually approaches the 
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FIG. 6: The T-fi phase diagram of square clusters near 
optimaUy doped (A*') ~ 3 regime ai U — 40 iUustrates the 
condensation of electron spin and onset of phase separation 
for spin degrees below Tp for various TV regions in fi space. 
The charge and spin susceptibility peaks, denoted by Tc and 
T* , define corresponding pseudogap regions calculated in the 
grand canonical ensemble, while boundaries for stable ferro- 
magnetic transitions fJ.+ {T) and fJ,-{T) are evaluated in the 
canonical ensemble. In equilibrium, the grand canonical value 
/ip at optimal doping (N) = 3 at T = reproduces the re- 
sult fj.p — {fi+ + fi-)/2 for canonical approach. The inset 
shows variation of canonical spin gap A" versus temperature 
for various values of /x. 

spin pseudogap as in the conventional BCS theory. 

Our exact calculations of phase diagrams in various bi- 
partite and nonbipartite clusters provide strong evidence 
for the existence of a narrow, homogeneous (pseudo)gap 
A* that vanishes near T^, coexisting with inhomoge- 
neous, weakly temperature dependent broad gap A"^, 
which disappears at higher temperatures > Tf. 
These phase diagrams display coherent and incoher- 
ent pairing (pseudo)gaps and possible superconductiv- 
ity in agreement with the recent STM measurements in 
HTSCs [l|,[3,ii,i,i,[3. 



B. Bipartite clusters at large U 

As for the cuprates, the bipartite clusters are useful 
for understanding the magnetic behavior and instabil- 
ities in manganites. The phase diagram in Fig. [5] for 
square clusters at J7 = 40, quite similar to other bipar- 
tite clusters at large U limit [2l[, displays characteris- 
tic features of managanites with strong electron corre- 
lations. In the ground state, the cluster at (N) — 3 
exhibits ferromagnetism in agreement with the Nagaoka 
theorem [36[. However, we observe saturated ferromag- 
netism, = I spin state for (N) « 3 clusters with one 



hole off half filling also at finite temperatures. The curve 
below Tp~^ signifies the onset of spontaneous magneti- 
zation with = I for parallel spin condensation. The 
positive charge gap (A'^ = 0.7787) for electron-hole (ex- 
citon) pairing manifests MH-like insulating behavior. In 
contrast, clusters show the minimum spin 5' = anti- 
ferromagnetism at (iV) w 2 and {N) « 4 in the ground 
state and at finite temperatures. The inset in Fig. [6] dis- 
plays the variation of spin gap for various regions. At 
/i = 1.35 the negative spin gap for {N) « 3 approaches 
zero as T ^ T^-^. Thus the region above T^~^ de- 
scribes a paramagnetic phase with zero spin gap for un- 
paired spins. In contrast, the positive spin gap in high 
doped regime at ^ = 0.25 changes its sign at tempera- 
tures above Tp- This picture describes a transition 
driven by temperature from antiferromagnetism into fer- 
romagnetism with S = 1. At half filling, a MH-like anti- 
ferromagnetism is stable at very low temperatures and 
unsaturated ferromagnetic state with S = 1 becomes 
more stable at higher temperatures, T > T^^^-^. How- 
ever, Xp'^^-'^ ^ as J7 — > oo and unsaturated ferromag- 
netism with S" = 1 at half filling can be stabilized at in- 
finitesimal temperatures. The well separated charge and 
spin susceptibility curves in the entire parameter range 
near {N) « 3 show spin-charge separation and decou- 
pling of charge and spin degrees. The susceptibility peak 
at T* for (A'^) « 3 in Fig. [S] displays a spin liquid behavior 
in the overdoped region for fj, < fj.p, while well developed 
negative spin gap in underdoped for ^ > /j,p regions at 
low temperatures describes a ferromagnetic insulator. In 
Fig. [HI the region of metallic-like behavior is manifested 
by the charge susceptibility peaks along the Tc curve. 

Phase diagram with fi dependent locally inhomoge- 
neous. A* < 0, and homogeneous. A" > 0, spin struc- 
tures at low temperatures, coexisting with charge or- 
dered homogeneous Mott-Hubbard like gap A'^ > dis- 
plays spin-charge separation and characteristic features 
of the CMR-manganite LaijjCajjMnOa and related mate- 
rials with alternating insulating ferromagnetic and charge 
ordered antiferromagnetic regions [El}. 



VI. CONCLUSION 

We have studied the dependence of the ground state 
and thermal properties in the repulsive Hubbard model 
on cluster size, geometry, electron number and interac- 
tion strength to understand the inhomogeneous super- 
conducting elements and stripes. The inhomogeneities 
found in exact calculations of clusters are promising for 
the description of geometric stripes with alternating su- 
perconducting and antiferromagnetic regions in high-Tc 
cuprates and magnetic domain structures in manganites. 
Spatial electron inhomogeneities capture the magnetic 
and pairing instabilities in clusters and respective bulk 
materials. The principal conclusion is that the exact 
solution for optimal inhomogeneities mimicked in small 
clusters can target essential features relevant to exist- 
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ing inhomogeneities in nanostructured materials on a 
nanoscale level. We found charge and spin gaps of equal 
amplitude in the ground state similar to the coherent 
pairing in conventional BCS theory. However, separate 
Bose condensation of electron charge and spin degrees 
with two consecutive transition temperatures into coher- 
ent pairing suggests a mechanism different from the pre- 
diction of the BCS coherent behavior with a unique crit- 
ical temperature. This picture is also consistent with 
the existence of two different energy scales for electron 
charge and spin pairing condensation temperatures in the 

HTscs flsiiai. 

The electronic instabilities in various geometries and 
in a wide range of U and temperatures will be useful for 
the prediction of coherent and incoherent electron pair- 
ings, ferroelectricity [1, and possible superconductiv- 
ity in nanoparticles, doped cuprates, etc. In contrast 
to bipartite clusters, the exact solution for the tetra- 
hedron depends on the sign of t and shows relatively 
weak dependence on U. For example, the tetrahedron 
exhibits similar features a,t U — 4 and 11 — 4:0 whenever 
t = 1. On other hand, the behavior of the tetrahedron 
for t = —1 strongly differs from that oi t = 1. These re- 
sults for frustrated clusters show that the properties are 
more sensitive to the change of the sign of the hopping 
term rather than the U parameter. This fact can explain 
why itinerant ferromagnetism can occur even at relatively 
weak interactions in frustrated systems. Our findings 
at small, moderate and large U carry a wealth of infor- 



mation regarding phase separation, ferromagnetism and 
Nagaoka instabilities in bipartite and frustrated nanos- 
tructures in manganites/CMR materials at finite temper- 
atures. These exact results allow us to understand the 
origin of level crossings, spin-charge sep aration, recon- 
ciliation and full Bose condensation m|. The obtained 
phase diagrams provide novel insight into electron con- 
densation, magnetism, ferroelectricity at finite tempera- 
tures and display a number of inhomogeneous, coherent 
and incoherent nanophases seen recently by STM and 
ARPES in numerous nanomaterials, assembled nanoclus- 
ters and ultra-cold fermionic atoms [l^, [s^ . 

Finally, we conclude that the use of the chemical po- 
tential and the departure from zero degree singularities 
in the canonical and grand canonical ensembles are essen- 
tial for understanding the important thermal properties 
and physics of phase separation instabilities and inho- 
mogeneities at nanoscale level. The article currently in 
progress is aimed to study the stability of pairing cor- 
relations and mag netism in the presence of transverse 
magnetic field [i^. It will be shown that magnetic flux 
tube inside the octahedral cluster can get trapped in sta- 
ble minima at half integral units of the flux quantum in 
hole-rich regions. 
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